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At low temperatures, dynamics in amorphous silicon occurs through a sequence of discrete activated 
events that locally reorganize the topological network. Using the activation-relaxation technique, a 
data base containing over 8000 such events is generated, and the events are analyzed with respect 
to their energy barrier and asymmetry, displacement and volume expansion/contraction. Special 
attention is paid to those events corresponding to diffusing coordination defects. The energetics is not 
clearly correlated with the displacement, nor with the defect density in well relaxed configurations. 
We find however some correlation with the local volume expansion: it tends to increase by about 4 
eV/A"^. The topological properties of these events are also studied; they show an unexpectedly rich 
diversity. 

PACS numbers: 61.20.Lc, 61.43. Dq, 66.30.-h, 67.40. Fd, 02.70.Lq 

Here, we apply it to an empirical model of amorphous 
silicon as described by a modified Stillinger- Weber po- 
tential The activation- relaxation technique allows 
one to concentrate on the activated mechanisms that are 
responsible for most of the dynamics below melting. We 
have already reported a first stage of this study in a recent 
Letter where we concentrated on only those mecha- 
nisms involving no coordination defects. Here we look at 
a wider spectrum of mechanisms, with a special emphasis 
on the diffusion of coordination defects. 

This paper is organized as follows. In section || we 
describe the activation-relaxation technique and give the 
details of the simulations that result in the data base 
of events. Next, in section H], we present the results 



I. INTRODUCTION 

The properties of amorphous semiconductors can vary 
widely as a function of the details of the preparation 
method: hot-wire, electron-beam deposition 0| and ion 
bombardment can yield samples with a significant 
spread in electronic and structural properties. Such di- 
versity is a reflection of the immensely large number of 
metastable configurations of nearby energy. 

The topological structure of this complex energy sur- 
face can be sampled indirectly, for instance by light il- 
lumination or ion bombardment, bringing a sample from 
one metastable state to another, often in a (statistically) 
reversible manner. A direct study of this energy surface 
requires the identification at the microscopic level of the 
mechanisms responsible for moving from one metastable 
state to another, and is much harder to perform. Because 
of the high degree of disorder, very few techniques can 
provide a truly microscopic representation of the bulk 
dynamics ||^. At best, one can extract some quantity av- 
eraged in time and space, providing a very rough picture 
of what is really happening. 

In spite of these difficulties, the past few years have 
seen significant experimental and theoretical efforts to 
try to provide some insight into the bulk dynamics. As 
it is becoming evident that little hard and precise in- 
formation about the local environment can be obtained 
by via static methods, more and more emphasis is put 
onto the development of techniques to sample dynamical 
quantities. 

This article presents the first detailed study of the mi- 
croscopic nature of activated mechanisms in a-Si. The 
method that we have used is the activation-relaxation 
technique (ART), introduced by us a few years ago [Q. 



extracted from this data base. Because this type of work 
is rather new, we also discuss the analysis of the data in 
some detail. The results contain both a global analysis 
and a more detailed topological classification. 



II. METHOD OF GENERATION OF THE EVENT 
DATABASE 

In this article, we concentrate on identifying and clas- 
sifying activated mechanisms that are responsible for the 
relaxation and diffusion in amorphous silicon. This is 
done using the activation-relaxation technique (ART), an 
energy-landscape method, that searches for barriers and 
new states in complex landscapes. 

As we shall see, the results we obtain are in general 
agreement with many of the experimental results men- 
tioned above and provide some bounds on the type of 
mechanisms that can take place in a-Si. Of course, em- 
pirical potentials have strong limitations, especially far 
away from the equilibrium position for which they are 
developed. Without giving too much weight to the exact 
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numerical values of the activation energies, it is neverthe- 
less possible to give a first broad picture of the wide vari- 
ety of mechanisms that can be associated with diffusion 
and structural relaxation. To go beyond the results pre- 
sented below, it will be generally necessary to use more 
accurate interactions such as tight-binding or plane wave 
methods. 

A. sampling one event with ART 

The aim of ART is to sample minimum-energy paths, 
starting in a local energy minimum, passing through a 
first-order saddle point (where the minimum-energy path 
has its highest point), and leading to another local en- 
ergy minimum. ART does this in three stages: leaving 
the so-called "harmonic well" , convergence to the saddle 
point, and convergence to the new minimum. The last 
stage, relaxation to a local energy minimum, is straight- 
forward and can be achieved by a wide range of standard 
minimization techniques, for instance the conjugate gra- 
dient (CG) method 0. The first two sta ges represent 
the activation to a saddle point and are specific to ART. 
As the end conditions of the first stage (leaving the har- 
monic well) are set by the actual implementation of the 
second stage (finding the saddle point), we will discuss 
these two stages in reversed order. 

The second stage is convergence to a first-order saddle 
point. At such a saddle point, the gradient of the en- 
ergy is by definition zero in all directions, and the second 
derivative of the energy is positive in all directions but 
one. The single direction with negative curvature is that 
along which the minimum-energy path proceeds. Within 
ART, we make the assumption that the direction of the 
minimum-energy path in the saddle point and the direc- 
tion towards the original local energy minimum have a 
significant overlap, i.e. that their dot-product is signif- 
icantly non-zero. If this assumption holds, a modified 
force field can be introduced in which the saddle point is 
a minimum This field is defined by the force G: 

G = F - [1 - a/(l + Ax)] (f • Ax) Ax, (1) 

where F is a SA^-dimensional force vector obtained from 
the first derivative of the potential energy. Ax is the 
displacement vector from the minimum, and a = 0.15 
a parameter determining how fast the motion is to the 
saddle point. In this second stage, the redefined force 
G is followed iteratively, usually along conjugate direc- 
tions, starting from just outside the harmonic well around 
the original local energy minimum. Ideally, this process 
would bring configuration directly to the saddle-point 
and stop there, but because the projection is an approx- 
imation of the valley, the configuration passes in the vic- 
cinity of saddle point without halting. We therefore stop 
the activation as soon as the component of the force F 



projected onto the displacement Ax changes sign, an in- 
dication that a saddle point has just been past. At that 
point, the activation is stopped, the configuration stored 
as the activated configuration and we move to the relax- 
ation. A more accurate convergence to the saddle point 
can be obtained by following directions determined by 
the eigenvectors of the dynamical Hessian |^ , but this is 
too costly for the system size of interest here. 

Most saddle points, even those belonging to ener- 
getically favorable minimum-energy paths, cannot be 
reached by following this redefined force G from a point 
well inside the harmonic region, i.e., the region of the 
energy landscape that is well approximated by a high- 
dimensional parabola centered on a local minimum. We 
refer the reader to Ref. jl^ for a detailed discussion of 
the origin of this problem and just describe our algo- 
rithm. We thus have to make sure that the configuration 
has left the harmonic well before following G. This first 
stage is implemented as follows. In a local energy mini- 
mum configuration, a few atoms and their nearby neigh- 
bors are displaced randomly. The total displacement is 
small, typically 0.01 A, and serves to create a non-zero 
force. At this point, the force will mostly point back 
to the minimum, and the component of the force, paral- 
lel to the displacement, will dominate the perpendicular 
components. However, we then increase the displacement 
from the original minimum until this no longer holds, and 
consequently we are outside the harmonic region. At that 
point, we can start the second stage of the activation and 
follow G. 

B. generating the data base 

The database of events that we use in the current 
manuscript is generated as follows. The energy land- 
scape is described by the Stillinger- Weber potential 
modified as described below. This empirical interaction 
includes a two-body and a three-body interaction: 

E^Y. ^(^«^- ) + E ^('^■'^- ' ' ^^-^fc) (2) 

where the brackets ( and ) indicate that each bond or 
angle is only counted once; the two-body potential is 

V{u,) = {Brif - 1) exp [{n, ~ a)-'] (3) 

and the three-body potential is 

V{nj,rik,Ojik) = eA (cos%fe + ^)^ 

X exp [7(ry - a)^^] exp [jirik - a)"^] (4) 

The numerical values for the parameters are A — 7.050, 
B = 0.6022, p = 4, a = 1.80, A = 31.5, 7 = 1.20, 
a = 2.O951A and e = 2.1682 eV; this set is identical 
to that used by Stillinger and Weber except for A which 
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has been increased by a factor of 1.5 in order to provide 
a more appropriate structure for a -Si |Q,|lO|;0 . Recent 
work on fracture underlines the fact that none of the 
available empirical interaction potentials describe silicon 
exactly |jl^ . This is particularly the case for energy barri- 
ers and density of defects; trends, more than exact values, 
are therefore what we are looking for here. 

We report here on results obtained from three indepen- 
dent runs. Each initial 1000-atom cell is constructed by a 
random packing in a large cubic cell. This configuration 
is then minimized at zero pressure to a nearby minimum 
state. ART is then applied iteratively with a Metropolis 
temperature of 0.25 eV in order to bring the configura- 
tion to a well-relaxed amorphous state; we consider that 
a system is "well-relaxed" when the energy does not de- 
crease significantly over hundreds of events. This takes 
place after about 5000 trial events, with a success rate 
slightly above 65%. 

After reaching a plateau in energy, for each event 
we store the initial minimum configuration, the saddle 
point configuration, and the final minimum configura- 
tion. Over the three independent runs, we collected a 
set of 8106 events. Figure |l| shows the radial distribu- 
tion function for run C at the beginning (CI) and end of 
the data acquisition run, 5000 trial-events later (C5000). 
Although radial distribution functions cannot discrimi- 
nate easily between realistic and non-realistic structure 
I p^ , the one obtained here is in good agreement with 
experimental data JI^ . 

Table | shows the structural properties of these two 
networks. A third of the bonds, involving 42% of the 
atoms, have been changed and yet the total energy and 
structural properties are almost unchanged. This is an 
indication that the configuration, through a sequence of 
events, has evolved considerably on an almost constant 
energy surface. 

TABLE I. Structural properties of configurations CI and 
C5000, the first and last structures used during the data col- 
lection in run C: the energy per atom in electron- Volt, the 
distribution of coordination using a cut-off at 3.05A, the av- 
erage bond angle and the width of the related distribution. 
The two configurations differ by 33% in their bonding and 
42% of the atoms have changed their neighbor list. 
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C5000 



Energy (eV) 
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e 

AO 



-3.903 
0.043 
0.948 
0.009 
3.966 

109.27 
10.04 



-3.918 
0.035 
0.960 
0.005 
3.97 
109.32 
9.74 



The bond angle distribution and coordination is com- 
parable with the best networks as obtained with realis- 
tic potentials |l^,|l^. A comparison with experiment is 
more difficult. Experimental works on amorphous silicon 
samples prepared by ion-bombardment report that ho- 
mogeneous samples of a-Si, without voids, have a density 
about 1.8 % lower than c-Si ||l|,|l3- Recent high-Q X- 
ray diffraction measurements on similarly prepared sam- 
ples show that well-relaxed a-Si could have an average 
coordination as low as 3.88, much below the 4.0 gener- 
ally considered to be appropriate for ideal o-Si . The 
authors of this work conclude that the high density of 
dangling bonds should be responsible for the lower den- 
sity of these samples. This analysis is not supported by 
our simulation, which leads to cells with a density as low 
as 7 % below that of the crystal, while keeping the num- 
ber of defects lower than that seen in this experiment. 
Clearly more work needs to be done to clarify this situ- 
ation, especially since the 12 % of dangling bonds seen 
by high-Q diffraction is at least an order of magnitude 
higher than what can be expected from either differential 
scanning calorimetry (1% defects) and electron-spin res- 
onance measurements (0.04% defects) These results 
also contrast with ab-initio and tight-binding computer 
simulations, which tend to find higher coordination, of- 
ten above 4.0 |l^]. The origin of this discrepancy is hard 
to identify at the moment but it can be due to a combi- 
nation of inaccurate interactions and/or the differences 
in time and lengths scales between simulations and ex- 
periments. In the case of our simulations, the stability of 
the configurations after 5000 trial-events suggests that we 
have reached some type of thermalization for this given 
interaction potential. 
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FIG. 1. Smoothed radial distribution function for config 
urations CI (solid line) and C5000 (dashed line) 
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III. GENERAL PROPERTIES OF EVENTS 
A. asymmetry, activation energy, and displacement 

The events can first be classified in term of three en- 
ergies: the energy asymmetry, the activation energy, and 
the total atomic displacement (see Fig. ||.) The sim- 
plest quantities to use for the classification of events are 
the barrier and asymmetry energies. In Ref. we give 
the distribution of these two quantities for the full set of 
8106 events. Both distributions are wide and relatively 
smooth. To push further the analysis, it is useful to es- 
tablish a first classification based on topological proper- 
ties of the network. 

The radial distribution function of relaxed a-Si goes to 
zero between the first and second neighbor. The middle 
of this region between first and second-neighbor lies at 
around 3.05 A. This allows us to establish a clear defi- 
nition of nearest-neighbor bonds between atoms. Struc- 
tural properties of the first and last minimum-energy con- 
figurations of run C in our database are listed in Table |; 
we find similar numbers for runs A and B. Most atoms 
that change neighbors in an event are four-fold coordi- 
nated both before and after. We discern three topolog- 
ical classes of events; if all atoms involved in an event 
keep their coordination unchanged between the initial 
and final state, we term it a perfect event; if topologi- 
cal defects change place during an event, but their total 
number is conserved, we have a conserved event, which 
describes defect diffusion; the remaining events are called 
creation/annihilation events. 

The exact details of the number of each class of events 
depends slightly on the cut-off radius between first and 
second neighbors; although low-energy configurations 
have a clear gap between first and second neighbor peaks 
in the radial distribution function, many saddle and some 
high-energy final configurations show structure in this 
gap. For consistency, we have chosen a fixed cut-off at 
3.05 A for all our analysis; the topological classification 
of some events might be affected by these parameters 
but the overall conclusions are not sensitive to the fine 
tweaking of that value. 

In a previous letter, we have discussed in some detail 
the class of perfect events (802 events) The bulk 
of our database comprises creation/annihilation events, 
with 5325 events, but until now we did not succeed in 
revealing interesting characteristics from these. In this 
paper, we focus on the 1979 conserved events in our 
database, describing directly the diffusion of defects with- 
out the creation or annihilation of coordination defects. 



Saddle 




Initial Final 

FIG. 2. The energy barrier and asymmetry in a schematic 
event as created by tlie activation-relaxation teclinique. 

As already mentioned, we produced 1979 conserved 
events, i.e., events where the number of coordination de- 
fects is identical in the initial and final state. The distri- 
bution of barriers and minimum to minimum energy dif- 
ferences is shown in Figure ^(a). The front of the barrier 
distribution peaks at about 4.5 eV while the asymmetry 
peaks at about 2.1 eV. This distribution is very similar 
to the total distribution shown in Fig. 1 of Ref. ||^. In 
fact, one of the most striking results we obtain is that the 
distribution of barriers and asymmetry is almost indepen- 
dent of the subset of events we select. The bottom box 
of Fig. ^, for example, compares the asymmetry distri- 
bution for the three classes of events: perfect, conserved, 
and creation/annihilation. Except for the peak at eV 
in the case of perfect events, involving atomic exchanges 
without modification of the overall topology of the net- 
work, the three distributions fall almost on top of each 
other at low asymmetry. The maximum of each distribu- 
tion is slightly shifted, however, with peaks at about 2.3 
eV for the conserved events, 2.5 eV for the perfect events 
and 3.2 eV for creation/annihilation events. 
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FIG. 3. Top: Distribution of the energy barrier (dashed 
hne) and energy assymetry (sohd hne) for the 1979 conserved 
events. Bottom: Comparison of the renormahzed asymmetry 
energy distribution for the conserved (dashed hne) and the 
perfect events (dotted hne) and all the others (solid line). 

The bias towards higher energies in the asymmetry 
distribution is expected: the distribution includes all at- 
tempted events, and not just those that are accepted; 
since the run is started in an already well-relaxed config- 
uration, most events lead to higher energy configurations. 

A similar insensitivity can be seen with the activation 
barriers. Although the precision on the barrier height is 
about 0.5 eV within the modified-force approximation, 
we can still say something about activation. Fig. ^ also 
plots the distribution of barriers. It peaks at about 4.0 
eV Taking into account the uncertainty on the empiri- 
cal potential, this result are consistent with experimental 
measurement. Shin and Atwater |2^] conclude, based on 
conductivity measurements, that the activation-energy 
spectrum extends from as low as 0.25 eV to about 2.8 
eV, supposing that a prefactor (entering logarithmically 
in the relation) is of order 1. Using isothermal calorime- 
try, Roorda et al find relaxation with a characteristic 
time of about 110 s between 200 and 500 C, also indi- 
cating a high activation barrier. Moreover, these results 
seem to depend only weakly on the method of prepara- 
tion (vacuum evaporation or ion implantation) and are 
limited by the fact that above 500 C the samples tend 
to crystallize. Both results are therefore also consistent 
with a continuous distribution of activation barriers. 

The tail of the distribution goes much beyond exper- 
imental values, and extends past 20 eV. Although such 
mechanisms are clearly unphysical, they underscore the 
fact that ART does not suffer from slowing down as the 
height of the barrier increases. This method is perfectly 
at ease with barriers of 0.1 eV as well as those of 25 eV. 
In the rest of this paper, we will concentrate on events 
with the 1 147 more physical barriers of less than 8 eV. 



Figure Q shows the histogram of the total displacement, 
defined as the square root of the sum of the square of each 
single- atom displacement, for these conserved events with 
a barrier of less than 8 eV. There is again little structure 
in the distribution. We note that the the average dis- 
placement to the saddle point is shorter than that to 
the new minimum. Based on preliminary simulations in 
other materials, this trend, although intuitively reason- 
able, is not always present and might be indicative is 
certain type of activation; more work remains to be done 
to clarify this issue. 
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FIG. 4. Distribution of total displacement at the barrier 
(dashed line) at at the new minimum (solid line) for conserved 
events. 



B. Volume expansion/contraction per event 

If is often suggested that there should be correlation 
between the energy of an event and its size. This im- 
mediately raises the question as to whether mechanisms 
by which the structure rearranges itself are local or non- 
local. 

Based on the observation that isothermal heat release 
curves obey bimolecular reaction kinetics and that the 
ion-beam-induced derelaxation scales with the density of 
displaced atoms due to the ion bombardment and ap- 
pears to be independent of electronic energy-loss mech- 
anisms, Roorda Q concludes that point-defect annihila- 
tion should control the structural relaxation. It remains 
unclear, however, whether this means the actual removal 
of defects or simply a clustering or a passivation by hy- 
drogen atoms. A similar relaxation (to within a factor 2) 
of ion-bombarded c-Si and a-Si suggests that both ma- 
terials have similar relaxation mechanisms but not nec- 
essarily identical. This similarity would point towards 
relatively local mechanisms of defect diffusion and relax- 
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ation since at this length scale, crystalline and amorphous 
materials resemble each other closely. 

Based on EPS data however, MuUer et al. suggest that 
from 1000 to 10 000 atoms have to move marginally in 
order to enable a single dangling-bond defect to anneal 
pl[ . This would be, at least qualitatively, in agreement 
with XPS measurements that suggests that the forma- 
tion of dangling bonds in a-Si:H under exposure to light 
is also accompanied by long-range structural rearrange- 
ments of the amorphous network [ p2[ , but it is in clear 
disagreement with the conclusion of Roorda et al. men- 
tioned above. The annealing mechanism suggested by 
this group, for example, is mutual annihilation of low- 
and high-density defects - vacancy/interstitial. Roorda 
et al. propose defects similar to that of c-Si but not 
necessarily identical Q . Mossbauer experiments with Sn 
suggest that vacancies exist also in a-Si [|3). Part of 
the difficulty in assessing more clearly the type of defects 
involved in relaxation and diffusion is that amorphous 
silicon crystallizes at about 500 C, rendering studies of 
self-diffusion very difficult |Q. 

Because it is not always clear what size means, we 
consider here three definitions: number of atoms, total 
displacement and local density deformations. 

The size of events is usually related to the number 
of atoms involved in the rearrangement of the network, 
i.e., the number of atoms that are displaced more than a 
threshold distance Tc. In Fig. || we plot the number of 
atoms involved as a function of rc, for the conserved and 
the full set of events. A local topological rearrangement 
will generally push the surrounding atoms outwards, or 
occasionally pull them slightly inwards. The distance 
Ar over which the surrounding atoms are pushed away 
(pulled inwards) will, because of elasticity arguments, 
scale with the distance r from the rearrangement as 
Ar r^'^ , for sufhciently large r. Alternatively, this 
can be rewritten as Ar ■ = Ve, where Ve is a constant 
volume independent of the distance r. The distance over 
which the atomic displacement Ar exceeds a threshold 
displacement rc is then equal to r = ^/Ve/Vc, and the 
number of atoms displaced more than Vc will then scale 
with Vr as gure shows that this 



scaling holds in the region where 0.1 < rc < 1.0 A and 
Ne ^ N. Selecting a threshold at 0.1 A the lower bound 
but also the typical vibration displacement at room tem- 
perature in Si, we find that, on average, about 50 atoms 
are involved in an event, clearly beyond the very local 
mechanism but well below the highest numbers proposed. 
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FIG. 5. Number of atoms displaced by a minimum thresh- 
old distance. The triangles and the X's are for saddle and new 
minimum position averaged over all events; the circles and 
diamonds are averaged over conserved events. In the case of 
activation, both distributions fall exactly on top of each other 
and it is hard to distinguish the X's from the diamonds. The 
solid line corresponds to the elasticity scaling, A'^ ~ r~^'^ as 
discussed in the text. 

The size of the events can thus be measured also with- 
out the introduction of a threshold distance, by measur- 



ing Ve = Ar 



averaged over some range of r. Figure 



shows the correlation between the event volume and the 
energy barrier and asymmetry. The figure shows that 
most proposed events tend to expand the sample locally 
by about Vg ^ 1 A^. Moreover, although scattered con- 
siderable, the data suggest a linear relation between the 
energy and the volume expansion with about 4 eV per 
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FIG. 6. Top; Asymmetry energy as a function of the vol- 
ume of the event (as defined in the text) for conserved events. 
Bottom: Distribution of the volume of events for conserved 
events at the saddle point (dotted line) and at the new mini- 
mum (solid line). 

We can also search for correlations between the total 
displacement and the asymmetry energy, which could re- 
late the diffusion length with the energy. This is plotted, 
for conserved events, in Fig. 0; there is little correlation. 
A similar negative result is obtained if we look at the 
correlation between the displacement and the activation 
energy. 
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FIG. 7. Asymmetry energy as a function of total displace- 
ment for the conserved events. 



events and the energy are difficult to establish due, in 
large part, to the wide spread of local environment typ- 
ical of disordered systems. These results can be used 
to put bounds on models of diffusion and relaxation in 
amorphous silicon. 



C. energetics of coordination defects 

Weak bonds and coordination defects form another re- 
curring theme in the study of dynamics and relaxation 
in a-Si. In this section, we discuss their properties in the 
sub-set of conserved events. The number of bonds bro- 
ken/created at the saddle point is 3.7 ± 1.4 and 3.4 ± 1.4 
with bond lengths of 2.43±0.09 and 2.62±0.12 A. At the 
final point, it is 4.3 ± 1.6 (equal number of bonds created 
and broken for a conserved event) with respective bond 
lengths of 2.46 ± 0.09 and 2.55 ± 0.10 A. 

These numbers are very similar to those obtained by 
concentrating on perfect events 1^. This reflects one of 
our main conclusions: that correlation between different 
properties of the events is weak. If the bond length of the 
resulting states is typically longer than that of the initial 
configuration, it is simply because the final configurations 
have typically a significantly higher energy. There is no 
preference for a single stretched bond; it is the medium- 
range strain that matters. 

One would expect a relation between relaxation and 
defect annihilation, but how strongly linked these two 
are is not clear a priori. Relatively little is know directly 
p5[ . To a first approximation, the total energy should 
follow the density of defects |2^ . Because of the similar- 
ity between crystalline and amorphous Si, Roorda et al 
concluded that relaxation occurs through defect annihi- 
lation 1^ . Polman et al find that Cu diffusion in annealed 
a-Si is 2 to 5 times faster than in non-annealed samples; 
this increase in diffusion might indicate a decrease in de- 
fects, trapping the Cu 

We see a certain correlation between defects and total 
energy as plotted in Figure ^(a). The distribution of en- 
ergy for a given number of defects is wide, going from 12 
eV for 35 defects to about 25 eV for 50 defects. Look- 
ing at the correlations of the energy with specific type of 
defects, we find much smaller impacts, with significant 
distribution in the total energy for a given number of 3- 
fold or 5-fold defects as is shown in the bottom panel for 
Figure |8[ 



The picture that emerges from these three approaches 
is that events are relatively localized, involving around 50 
atoms, and require some local expansion to take place, 
as would be expected. Correlations between the size of 
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FIG. 8. Top: Total energy as a function of the number 
of 4-fold atoms; bottom: total energy as a function of the 
number of 3-fold (open circle) and 5-fold (full diamonds) for 
all events. 

Clearly, therefore, the definition of defects must also in- 
clude strained environment and not just the coordination 
defects mentioned here — relaxing some highly strained 
ring might require the creation of a bond defect. Once 
again, the situation is much less clear than is generally 
thought. We must emphasize here that for less relaxed 
samples, the correlation between energy and the number 
of coordination defects is much better. 



D. Topological classification 

To go beyond the scalar picture given above, we need 
to consider in more detail the nature of the topological 
changes. The classification scheme applied here is an 
extension to the classification scheme that we used for 
the perfect events in an earlier letter All atoms that 
change their neighbors are alphabetically labeled, the 
topological change is determined by specifying the list of 
all bonds before the event, and of all bonds after. 



For perfect events, there are as many bonds before as 
after the event, and moreover, the set of all these bonds 
can always be organized into a ring of alternating created 
and destroyed bonds. This ring can be represented by the 
sequence of atoms visited. For instance, in event abacbd, 
the bonds before and after the event are ab, ac, hd, resp. 
ba, cb, da. Thus, bonds ac and bd are replaced by ad and 
be, while bond ah is present both before and after the 
event. Many equivalent rings exist, but the convention 
of always using the alphabetically lowest label makes this 
classification unique. 

If the event is not perfect, this classification scheme has 
to be modified. In case a single bond (or dangling bond) 
jumps from one atom to another, the set of bonds does 
not form a closed ring, but an open chain of alternating 
bonds before and after the event; the same classification 
scheme can still be used, with the note that the first 
and last atom are not bonded. It also happens that the 
event comprises a series of bond exchanges which are in 
disconnected regions (typically still nearby - interacting 
via the strain field). For such events, we introduce the 
concept of "ghost bonds" (represented by a dot in the 
label), that are added to either the set of bonds before 
or after the event. 

Using this topological analysis, we can have a first 
crack at the events. 1148 events in our database have 
an activation barrier of less than 8 eV, the other events 
might be considered unphysical. Of these, 447 are too 
complicated (involving too many defects or too many 
disconnected regions) to be analyzed, leaving 701 labeled 
events. 

For conserved events, there are no dominant labels, 
contrary to what is found for perfect events where three 
labels account for 85% of the events. Such diversity un- 
derscores the difficulty in trying to identify mechanisms 
and relating them to experimental information. Clearly, 
the dynamics of defects in amorphous silicon is much 
more complicated than is usually thought. 

In the large set of labels, an often occurring theme is a 
ring of bonds corresponding to the Wootcn- Winer- Weaire 
(WWW) bond-exchange mechanism H^,^ , which in our 
classification scheme has the label abacbd. We find that 
up to two WWW moves can take place in a single event. 
This rearrangement changes the local ring structure and 
redistributes the strain, affecting the jump barrier seen 
by the dangling bond. 

After removal of all these rings from the events as such, 
what remains is often a single coordination defect that 
jumps to a 1st, 2nd or higher neighbor. Table |l| presents 
a partial list of such labels. It is remarkable that the 
diffusion of coordination defects requires, in general, a 
topological rearrangement more complex than one would 
expect from the displacement of the bonds. Only nine 
events are of the abc type, the smallest rearrangement 
possible for the motion of a bond. 

Other events involve longer jumps, at least in topolog- 
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ical terms. The abacbde, for example, reflect this type of 
behavior. We show one such event in Fig. |[ Very few 
of the classified events displace more than one defect. (It 
could be that in the non-classified events this situation 
occurs more often.) 

IV. CONCLUSIONS 

Based on an extensive list of events in a-Si, represent- 
ing a wide range of characters, we can provide a general 
overview of the nature of the microscopic mechanisms 
responsible for the diffusion of topological defects in this 
material. To do so, we havei concentrated on a class of 

TABLE IL Most common events for low energy conserved 
events. The first column gives the classification label in terms 
described above. The second and the third columns refer to 
the number of events carrying this label and the number of 
topological defects — typically dangling bonds — displaced 
in the event, respectively. The last column, finally, gives a 
topological indication of how far the defect diffuses. A lo- 
cal jump brings a defect from one atom to its near-neighbor; 
non-local events will involve jumps to the second-, third-, or 
even fourth-neighbor. 
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local/non-local 


abacadefbfgbedc 


2 


1 


local 




abacadefegfbedc 


2 


1 


local 




abacb 


9 


1 


local 




abacbdaeb 


5 


1 


local 




abacbdaebfagb 


6 


1 


local 




abacbdaebfagc 


2 


1 


local 




abacbdaebfbgf 


4 


1 


2nd nei 


jhbor 


abacbdaec 


8 


1 


local 




abacbdaef 


6 


1 


local 




abacbdaefageg 


2 


1 


local 




abacbdbed 


5 


1 


2nd nei 


jhbor 


abacbdbedcfcfg 


6 


1 


local 




abacbdbedcfcgdf 


2 


1 


local 




abacbdbefegfdg 


4 


1 


local 




abacbde 


3 


1 


2nd nei 


jhbor 


abacdcedb 


8 


1 


local 




abacbdedfeb 


5 


1 


local 




abacdcedbfdfe 


2 


1 


2rd neif 


;hbor 


abacdcedbfdfg 


8 


1 


3rd neif 


;hbor 


abacdcedbfdfgeg 


4 


1 


3rd neif 


;hbor 


abacdcedfcb 


4 


1 


local 




abacdcefbfgbechdi 


3 


1 


3rd neif 


;hbor 


abacdebefbdcg 


11 


1 


2nd nei 


jhbor 


abacdefegfdcb 


3 


1 


local 




abc 


9 


1 


local 




abcde 


2 


1 


4th neif 


;hbor 


abcbdcebf 


7 


1 


local 




abcbdcebfcg 


5 


1 


3rd neif 


;hbor 


abcdc.ec.fdg.ehi 


1 


4 


2nd nei 


ghbor 




FIG. 9. Event of type "abcbdce". Barrier: 2.8 eV, Asym- 
metry: 2.1 eV, total displacement 2.5 A, four bonds broken 
and four created. 

events that involve the displacement of coordination de- 
fects while keeping their overall number constant. 

Analyzing a wide range of structural and topological 
properties of these events we find that: (1) In a well- 
relaxed sample, there is little correlation between the 
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number of defects and the total energy; the relaxation of 
strain can take place in more subtle ways, sometimes in- 
volving the creation of topological defects. (2) Taken into 
account the use of an empirical potential, the activation 
barriers are in agreement with experimental value. (3) 
We find little correlation between the activation barrier 
or the asymmetry and the deformation of the network, 
either in terms of the number of atoms involved or the to- 
tal displacement; the energy is best described in term of 
the volume of an event. (4) A topological analysis of the 
conserved events show an unexpected richness; we find 
literally hundreds of different mechanisms that cannot 
easily be put in a few classes. As a rule, the Wooten- 
Weaire- Winer bond exchange mechanism, dominant for 
perfect events, still plays a major role. Defect diffusion is 
often local, with coordination defects jumping from one 
atom to a neighbor, but it can also go as far the fourth 
neighbor, in a chain-like fashion. 

These results can be used to put bounds on models of 
diffusion and relaxation in amorphous silicon. For exam- 
ple, the Fedders and Branz model for a-Si:H states 
that (1) relaxing the defect structure often requires sev- 
eral atoms to move simultaneously, (2) only by cooper- 
ative motion do the position changes lower or conserve 
the total energy, (3) the size of the barrier generally in- 
creases with the number of atoms that must move simul- 
taneously. Our results support points (1) and (2) but not 
(3). ^ 

This study represents only a first step in the study of 
microscopic activated mechanisms in a-Si. More work 
remains to be done to converge barriers using more ac- 
curate interaction potentials. It is important also to try 
to connect some of these results with hard experimental 
numbers, a challenge both for theorists and experimen- 
talists. Already, however, we can see that the dynamics 
of disordered materials promises to be much more com- 
plicated than was thought before. 
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